Origin of n- and p-type conductivity in undoped «-PbO: role of defects 
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The first principles calculations (GGA) have been applied to study the crystallographic defects in 
a-PbO in order to understand an origin of n- and p-type conductivity in otherwise undoped a-PbO. 
It was found that deposition in the oxygen-deficient environment to be defined in our simulations 
by the Pb-rich/O-poor limit stimulates a formation of the O vacancies and the Pb interstitials both 
to be characterized by quite low formation energies ~ 1.0 eV. The O vacancy being occupied by two 
electrons shifts a balance of electrons and holes between these two defects to excess of electrons (four 
electrons against two holes) that causes the n-type doping. For the Pb-poor/O-rich limit, an excess 
of oxygen triggers a formation of the O interstitials characterized by such a low formation energy 
that spontaneous appearance of this defect is predicted. It is shown that the concentration of the 
O interstitials is able to reach the extreme magnitude equal to number of the possible defect sites 
(~ lO'^^cm"^). The localized state formed by the O interstitial is occupied by two holes and because 
there are no other defects in reasonable concentration to balance the hole redundancy, p-type doping 
is induced. 
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Recent developments in radiation medical imaging 
technology imposes new requirements on the photocon- 
ductive m.aterials: polycrystalline and amorphous mate- 
rials are considered for application in the flat panel de- 
tectors [l[ with preference given to high Z (atomic num- 
ber) polycrystalline semiconductors. Among all poly- 
crystalline materials, the Lead Oxide (a-PbO) is the 
most promising candidate due to its long history of em- 
ployment in optical imaging devices |31- However, de- 
spite such extensive application, properties of a-PbO are 
not well established but the low mobility-time product 
Q indicates the poor carrier transport 3]. Therefore, it 
is believed that the carrier transport in PbO compound 
should be controlled by trapping on defects. For appli- 
cation of PbO in real time procedures as fluoroscopy, the 
transport characteristics are required to be optimized for 
which an understanding of defect nature is crucial. The 
thermal evaporation technique used for deposition of the 
PbO layers [l|,|^ presumes a deficiency of the atomic oxy- 
gen ultimately leading to formation of the oxygen vacan- 
cies. The conductivity enhancement of the PbO samples 
observed after the thermal annealing in atmosphere of 
oxygen supports this hypothesis 4] . In our previous stud- 
ies devoted to defect formation in a-PbO single crystal 
we confirmed that this compound tends to loose oxygen 
due to the low formation energy of the oxygen vacancy 
0. The oxygen vacancies are found to be occupied by 
two electrons that led to a prediction of excess of elec- 
trons. The observed experimentally n-type conductivity 
of the PbO layers grown under conditions of the oxygen 
deficiency 0, m, M, 01 verifies our theoretical findings ^] . 

The samples being characterized by the oxygen excess 
show the exact opposite behavior: p-type conductivity 
6|. The Pb vacancies considered in our previous study 
5| , although they are occupied by two holes and the oxy- 
gen reach conditions promote their formation, can not 
be responsible for p-type doping because of their rather 
high formation energy. In order to develop the post- 
growth annealing procedures directed on suppression of 



the vacancy defects, an understanding whether origin of 
the p-type conductivity is being defect related is urgently 
required: application of the annealing in atmosphere of 
oxygen become questionable provided it generates new 
defect while healing the oxygen vacancies. Therefore, to 
receive a real picture the comprehensive study of all type 
of defects to be formed under the different growth condi- 
tions are required. Here we applied the first-principles 
methods for such study (the Pb-rich/O-poor and Pb- 
poor/O-rich limits to the growth conditions are consid- 
ered) with focus on an analysis of the defects being in- 
herent to the a-PbO crystal structure, their appearance 
inside the band gap, the charged states and the formation 
energies. 



I. METHODS 

We study a formation of the defects in a-PbO with 
help of the density functional theory (DFT) provided 
by WIEN2k package j8|] which adopts the full-potential 
augmented plane-wave -(- local orbitals method. The 
generalized gradient approximation (GGA) [9] with the 
Perdew-Burke-Ernzerhof parametrization was used. The 
5p, 5d, 6s and 6p electrons of the Pb atom and 2s and 
2p electrons of the O atom have been treated as the va- 
lence electrons (lower energy electrons were assigned to 
the core states). We applied the supercell approach of 
108-atom size (3x3x3 array of the primitive unit cells). 
The geometry optimization procedure applied to super- 
cell containing a defect was performed based on mini- 
mization of the forces [1(J| {RKmax=^)- The Brillouin 
zone was covered by 5x5x4 k-point mesh utilizing the 
Monkhorst-Pack scheme. 

The tetragonal lead oxide a-PbO is characterized by 
the layered structure to appear in polycrystalline form 
upon compound growth 1]. GGA tends to overestimate 
the interlayer separation in the layered structures while 
underestimating the band gap size. For a-PbO system 



in which the orbital overlap between layers controls a 
size of the band gap, it results in compensation effect !5|. 
In particular, for the lattice parameters optimized with 
GGA the band gap is only slightly underestimated as 1.8 
eV against the experimental value of 1.9 eV [11] (appli- 
cation of the experimental value of the lattice parameters 
causes the band gap to shrink by 0.22 eV [5]). Because 
a size of the band gap is crucial in order to define cor- 
rectly an appearance of the defect state inside the band 
gap, the band structure properties have been evaluated 
for the optimized lattice parameters. However, for calcu- 
lations of the defect formation energy, an adjustment of 
the interlayer distance by using the experimental data of 
the lattice constant cq = 5.07 A |l2| has been performed. 
It was found that only formation energy of the Pb in- 
terstitial is considerably affected by mismatch in lattice 
parameters to be caused by poor incorporation of Pb 
atom being of the large atomic radius between the layers 
(the formation energies of the vacancy defects have been 
indicated not to be affected by GGA limitation [5|). 
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II. CRYSTALLOGRAPHIC DEFECTS IN a-PBO 

Analysis of literature on defect formation in semicon- 
ductors [iJl and particularly in simple oxides [ij, has 
revealed that additionally to the vacancy defects inves- 
tigated in our previous study [5], the interstitials have 
to be considered. For the crystalline materials, the in- 
terstitial defects induce the large lattice distortion into 
the immediate neig hbourhood that results in its high for- 
mation energy [l3|. In this regard the crystal lattice of 
the polycrystalline a-PbO is unique, because distance be- 
tween weakly interacting layers is large enough to accom- 
modate the foreign atom (cq = 5.07 AfT^). Therefore, 
along with the vacancy defects formation of the intersti- 
tials is expected to be inherent to the a-PbO lattice. 

We start from short description of the vacancy defects 
being investigated in our previous work [5!]. The O va- 
cancy {V'~') is occupied by two electrons localized on the 
four nearest Pb atoms with the electron density maxi- 
mum to manifest at the vacancy site (see the electron 
density distribution in Fig. [T] (a)). The O vacancy can 
appear in three charged states (0, 1+, 2+). The local- 
ized state generated by the O vacancy is located deep 
inside the band gap at Ed{V'^)-Ev=1-03 eV, where Ey 
is the top of the valence band. In contrast, the Pb va- 
cancy creates the localized state occupied by two holes 
and the defect level appears close to the valence band top 
at EDiV^^)-Ev=0-l eV. The electron density is spread 
between O and Pb atoms showing the long defect tails 
(see Fig.[T](b)) that is not surprising taking into account 
the defect level appearance close to the valence band top 
(see Fig. [2]) . The Pb vacancy has also three charge states 
(0, 1—, 2—) but its charge states would be of negative 
sign, i.e. ionization of the vacancy occurs through elec- 
tron gaining. 

There are two interstitial defects to be considered in a- 



FIG. 1: Schematic presentation of the crystallographic de- 
fects in the tetragonal Q-PbO (space group 129P4/nmm) and 
the electron density map plotted with XcrySDen [la ] for the 
energies Ed ±0.05eV^. (a) the O vacancy, (b) the Pb vacancy, 
(c) the O interstitial and (d) the Pb interstitial. 
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FIG. 2: Energetic location of the states induced by the vari- 
ous defects (1/°, V^"" , 7°, 7^'^) within the band gap of Q-PbO 
(_Bg=1-8 eV). Electrons occupying the vacancy site are indi- 
cated by the arrows while filled circles correspond to absence 
of electron. 



PbO crystal lattice, the O and Pb interstitials. Binding 
of the O interstitial (/'"') to the host occurs to single Pb 
atom by making use of its lone pair PbtGs" electrons. The 
lowest energy state is achieved when the O interstitial 
is positioned just above the Pb atom (see Fig. [TJc), the 
bond length is 2.03 A). Inside the band gap, the localized 
state is located above the midgap at Eu{l'~^)-Ev—l-27 
eV (see Fig. [J). This localized state can accommodate 



up to two extra electrons and, therefore, there are three 
charge states (0, 1 — , 2—) are assigned to it. 

The energeticahy favourable position of the Pb inter- 
stitial (/^^) in the lattice is similar to the O intersti- 
tial, i.e. on top of the Pb atom. The Pb interstitial 
is exited from Pb:6s^ 6p^ state to the Pb:6s^ 6p^ state. 
The host Pb atom participates in bonding again with its 
Pb:6s^ electrons (the bond length is 2.9 A). The origi- 
nated defect state is occupied by two electrons which are 
Pb:6p^ electrons belonging to I^^ (see Fig.[T](d)). Since 
each electron stays on its own orbital, the Hund's rule 
dictates the ferromagnetic ordering of these electrons. 
The Pb interstitial induces the localized state inside the 
band gap and spin-polarized solution separates the occu- 
pied and unoccupied states as E£){I^^)-Ev=0.99 eV and 
Ed{I^^)-Ev = 1-51 eV (see Fig. [21), respectivelly. Ioniza- 
tion of the defect state formed by the Pb interstitial can 
occur as through electron gaining ((1 — ) for one electron 
added and (2—) for two electrons) so through an excita- 
tion of existing electrons. Overall, there are five charged 
states (2-I-, 1-I-, 0, 1 — , 2—) can be considered for the Pb 
interstitial. 

Therefore, we have established that some defects can 
act as a donor impurity (V^^o) -^ yO(i+/2+)^^ q^j^_ 
ers as acceptors {V^^^°^ -^ yPb(i-/2-) ^^^ jO(o) ^ 
/0(i-/2-)-|^ while I^^ can be both an electron donor 
(jPb(o) ^ ^Pb(i+/2+)) ^^^ ^^ acceptor (jPb(o) ^ 
/Pb(i-/2-)) An ability of defects to appear in the several 
charge states can promote the electron exchange between 
defects: the electron donor compensates missing electron 
on the acceptor. The process of the charge exchange de- 
fines the energetically favorable charge state for defects 
to appear upon compound growth. 



III. FORMATION ENERGY OF DEFECTS 

The main way to distinguish a formation of which de- 
fect is inherent to the crystal structure is to compare 
the formation energies of different defects as particularly 
this parameter is a measure of the defect concentration. 
The formation energy of a defect D in charge state q 
(-h2/-hl/0/-l/-2) can be defined as [Hi 

AEf{D) = EtotiD''yEtotiSyJ2'^'f''+1^^F+^'''+^^^ 

i 

(1) 

where Etot{E>'^) and Etot{S) are the total energy of the 
system containing the single defect and defect-free sys- 
tem, respectively, rii indicates a number of i-atoms re- 
moved (rii > 0) or added (rii < 0), while ^i is the 
chemical potentials. {Ep + Ey + AV^) is the position 
of the Fermi level relative to the valence band maximum 
{Ev) which has to be corrected by /S.V (for details see 
Ref. [13]). Another correction term, so-called band gap 
error A.Eq has been applied as well [ig. To simulate 
the deficit and excess of oxygen, the Pb-rich/0-poor and 
Pb-poor/0-rich limits have been applied in evaluation 



(a) Pb-rich / 0-poor (b) Pb-poor / 0-rich 



7.0- 


(0/1 -) 
V^-/2-) : 


> 6.0- 

^ 5.0- 
ho 


NsV'' 


u 4.0- 


I ; (0/2-) ^ 


c= 3.0 

o 

■■§ 2.0- 

— < 


I^^ (0/1 -) Nv 


o 1.0- 




^p^'irmK 


0.0- 
.1 n- 


v° 







:(0/1-) 


2-) 


4.0- 


(0/1 -) 

yi-/2^ 


. \^" 




^(im\ 


,v° 


3.0- 


As 




\ 




(2+/1+)^ 




tPBS 


2.0- 








1.0- 






^Pb 


00- 


1° 


'' ini9-\ 


\ 


1.0- 




K 



0.0 0.6 1.2 1.8 0.0 0.6 1.2 1.8 



Ep,eV 



Ep,eV 



FIG. 3: Formation energy as a function of the Fermi level posi- 
tion {Ef varies from to -Eg^I.S eV [l3)- (a) Oxygen deficit 
(Pb-rich/O-poor Hmit). (b) Oxygen excess (Pb-poor/0-rich 
limit). 



of the chemical potentials. Using A/iJ°(PbO)=-2.92 eV 
[a|, the chemical potentials have been found for the Pb- 
rich/O-poor limit as //(o)=-Eto*(0 atom)-|-A/iJ°(PbO) 
and ^(Pb)=£'tot(Pb atom), while for the Pb-poor/O- 
rich limit as /i(o)=-Etot(0 atom) and /i(pb)=£'tot(Pb 
atom)+A^7J°(PbO). 

We have calculated the formation energies of the differ- 
ent defects in the relevant charge states and our results as 
a function of the Fermi level position for the Pb-rich/O- 
poor and Pb-poor/0-rich limits are shown in Fig. [3] As 
Hi deviates in a range defined by A/i?°(PbO) (in real 
growth conditions to be controlled by the partial pres- 
sure of reactants) it induces a shift of the formation ener- 
gies. The concentration of the defects can be estimated 
through their formation energies as Ndefects ^ NsUes 
exp{—AE^{D)/kbT), where ki, is the Boltzmann con- 
stant, T — 570 K is the deposition temperature and 
N sites is the concentration of sites available for defects 
to be formed. 

For the Pb-rich/O-poor limit (see Fig. |3](a)), the for- 
mation energies of the O vacancy (I^*"*) and the Pb inter- 
stitial (/^'^) are appeared to be much lower than that for 
other defects and, therefore, these defects are expected 
to dominate. The defect concentrations estimated for 
their zero charge states are lO^^cm^"^ for the O vacancy 
and 10^^cm~'^ for the Pb interstitial. For the Fermi level 
assigned to the midgap, the (0, 1-|-) and (0, 1-) states 
are energetically preferable for the O vacancy and the 
Pb interstitial, respectively. Since the thermodynamic 
transition y'-'(o) -^ yO(i+) is found to occur almost at 
the midgap, it indicates a potential for both states to be 
formed during the compound growth. However, to be 



ionized, y'-'^") has to dispose its electron and therefore, 
appearance of the T/0(i+) charge state would naturally 
depend on presence of the electron acceptor. For the Pb 
interstitial, we have calculated the formation energy for 
the five states (2+, 1+, 0, 1 — , 2—). Since the forma- 
tion energy of (1-) state to occur when an extra electron 
is added is significantly lower than that when electron 
is removed {AE^ {I^^'^^+^)=3.98 eV for Ep at midgap), 
the (1-I-) and (2+) states have been disregarded. The 
thermodynamic transition J^^C) — >■ /Pt>(i-) occurs close 
to the midgap and, therefore the charge states (0, 1—) 
appear to be the preferable states. To become ionized 
to the (1—) state, /P'^(o) has to gain electron (electron 
acceptor). 

Therefore, the charge exchange between T/'-'(o) and 
/Pb(o) jg possible as one of them can work as an elec- 
tron donor while other as an acceptor. The concentra- 
tion of the ionized states manifesting upon the charge 
exchange would be limited by defect of lower concentra- 
tion, i.e. by the Pb interstitial (lO^^cm"'^). It is impor- 
tant to emphasize, that a balance of electrons and holes 
occupying both defect sites in either combination of the 
charge states (i/0(o)/jPb(o) ^^ v^o(i+)/jPb(i-)) jg always 

shifted to a higher population of electrons (four electrons 
against two holes). Therefore, in respect to a conductiv- 
ity type, redundant electrons induced by the O vacancies 
would cause n-type doping of a-PbO at the Pb-rich/0- 
poor limit to be in agreement with our previous studies 
[5|. In the non-equilibrium electron gas, the ionized O 
vacancies T/*-'(i+) and the Pb interstitials in both charge 
states /P'^f"), jPt'Ci-) would work as the trapping centers 
for electrons in the conduction band. 

For the Pb-poor/0-rich limit, there is only the O inter- 
stitial, /°, to possess the low formation energy, while for- 
mation energies of other defects {I^^, V'~^ and V^^) are 
higher than 3.0 eV. In fact, the formation energy of /*-* 
being so closed to zero implies the spontaneous appear- 
ance of this defect such as the defect concentration can 
become equal to the concentration of the possible defect 
sites {Nsites ^ 10^^cm~^). The O interstitial in its neu- 
tral charge state /'-'C) misses two electrons and to be ion- 
ized to the charge states (1-/2-) has to acquire electrons. 
Since the thermodynamic transition /'-'(o) -^ /0(2-) oc- 
curs close to the midgap, l'^ potentially can be formed 
in the charge states (l-,2-). However, because concentra- 
tion of centers of the donor type {I^^, V^) is negligible, 
the O interstitial is expected to occur primary in its l'~'^^> 
state. The /°(°) state would act as a trap for the non- 
equilibrium electrons in the conduction band. Therefore, 
under growth conditions corresponding to the oxygen ex- 
cess, the high concentration of the /'-'('') defects each to 
be occupied by two holes, would be able to cause p-type 
doping. Moreover, the state formed by the O interstitial 
is characterized by quite delocalized nature (see Fig. [1] 
(c)): the long defect tails in combination with the high 
defect concentration would results in strong defect-defect 
interaction potentially leading to formation of the defect 
band. 



It is also important to emphasize that all investigated 
defects, V'^ , V^^, l'^ and I^^ being considered in bulk 
a-PbO would be inherent to surface of the single crys- 
tal (platelet). For some defects, the formation energy is 
found to be lowered when they are placed on the platelet 
surface. The strongest effect is observed for the Pb in- 
terstitial for which the formation energy drops by '^ 1.0 
eV. For the vacancies, the formation energy is only af- 
fected provided they appear at the edge of the platelet 
(for example, it is reduced by 0.3 eV for oxygen vacancy). 
The drop in the formation energy at the platelet edge is 
a result of the surface reconstruction leading to forma- 
tion of the Pb-0 double bonds instead of single bonds. 
Since a balance of electrons and holes occupying the de- 
fect sites is not affected by formation of the surface states 
then considered here phenomenon of the n- and p-doping 
remains to be valid. There is one more interesting fea- 
ture to be examined which is a defect pairing. For the 
Pb-rich/0-poor conditions, the pairing may occur be- 
tween V'^ and I^^ defects to be promoted by lowering 
in the formation energy by 0.84 eV. Such pair V''^-!^^ 
still forms two defect levels at Ed{V^^)-Ev^1-22 eV and 
Ed{V^^)-Ev=0.66 eV and each level is occupied by two 
electrons. For the Pb-poor/0-rich conditions, the inter- 
action between 1° and V^^ potentially can occur but be- 
cause each localized state is occupied by two holes, such 
interaction is repulsive. Under the equilibrium growth 
conditions, the formation energy of V^^ and V would 
become comparable indicating a possibility of their pair- 
ing as well tl8|. This combination of defects is unique 
as it induces vanishing of the dangling bonds and, there- 
fore, does not generate the defect state inside the band 
gap. Moreover, because the formation energy is lowered 
by 1.47 eV upon defect pairing, it indicates the pairing 
to occur not only during the compound growth but also 
as a post-growth migration of the O vacancy towards the 
Pb vacancy. 



IV. CONCLUSION 

Our study on formation of the defects in single crys- 
tal a-PbO has revealed that some defects under specific 
growth conditions would appear in significant concentra- 
tion causing the n- or p-type doping of compound. In 
particular, for the Pb-rich/0-poor limit a formation of 
the V'^ and I^^ defects are found to dominate with con- 
centrations of lO^^cm"'^ and lO^^cm""^, respectively. The 
obtained defect concentrations are in good agreement 
with concentration of the localized states estimated ex- 
perimentally (10^'*-10^'^cm~^ 2]). It was estabhshed that 
accumulation of electrons on sites of the V'~' and I^^ de- 
fects results in n-type doping of a-PbO. For the oxygen 
excess created during the compound growth (Pb-poor/0- 
rich limit), another defect is found to dominate which is 
the O interstitial. The concentration of /^ is predicted 
to approach the highest limit because its formation en- 
ergy is around zero, i.e. defect concentration may be- 



come equal to number of the possible defect sites. Since 
each j'-'(°) is occupied by two holes, its high concentra- 
tion would result in p-type doping of a-PbO. Therefore, 
because concentrations of l'^ , V'^ and I^^ defects induc- 
ing the defect levels close to the midgap are considerably 
high, their impact on the transport properties of a-PbO 
would become apparent. Although, the post-growth an- 
nealing is an efScient procedure which can be applied to 
reduce concentration at least of the vacancy defects, but 
there is always a high risk to worsen the transport prop- 
erties of PbO through generation of the new defects such 



as intcrstitials. 
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